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ABSTRACT 

We discuss the second order contributions to lensing statistics resulting from the clus- 
tering of background sources from which galaxy shape measurements are made in weak 
lensing experiments. In addition to a previously discussed contribution to the lensing 
skewness, background source clustering also contributes to the two-point correlation 
function, such as the angular power spectrum of convergence or shear. At arcminutc 
scales or above, the second order contribution to the angular power spectrum of con- 
vergence due to source clustering is below the level of a few percent. The background 
clustering of sources also results in a non-Gaussian contribution to the power spec- 
trum covariance of weak lensing convergence through a four-point correlation function 
or a trispectrum in Fourier space. The increase in variance is, at most, a few percent 
relative to the Gaussian contribution while the band powers are also correlated at the 
few percent level. The non-Gaussian contributions due to background source cluster- 
ing is at least an order of magnitude smaller than those resulting from non-Gaussian 
aspects of the large scale structure due to the non-linear evolution of gravitational 
perturbations. We suggest that the background source clustering is unlikely to affect 
the precision measurements of cosmology from upcoming weak lensing surveys. 
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1 INTRODUCTION 

Weak gravitational lensing of faint galaxies probes the dis- 
tribution of matter along the line of sight. Lensing by large- 
scale structure (LSS) induces correlation in the galaxy ellip- 
ticities at the few percent level and can be detected through 
challenging statistical studies of galaxy shapes in wide-field 
surveys (e.g., Blandford et al. 1991; Miralda-Escude 1991; 
Kaiser 1992). An important aspect of weak lensing is that 
these ellipticity correlations, or associated statistics, provide 
cosmological information that is, in some cases, complemen- 
tary to those supplied by the cosmic microwave background 
data at the same level of precision or better (e.g., Jain & Sel- 
jak 1997; Bernardeau et al. 1997; Kaiser 1998; Schneider et 
al. 1998; Hu & Tegmark 1999; Cooray 1999; van Waerbeke 
1999; see Mellier 1999 and Bartelmann & Schneider 2000 
for recent reviews). Indeed, a wide number of recent stud- 
ies have provided the clear evidence for weak lensing due to 
large scale structure (van Waerkebe et al. 2000; Bacon et 
al. 2000; Wittman et al. 2000; Kaiser et al. 2000), though 
more work is clearly needed to understand even the statis- 
tical errors and biases. While biases can result from certain 



aspects related to the selection of observational fields (e.g., 
Cooray et al. 2000), statistical errors include those that are 
fundementally present due to the non-Gaussian nature of 
the large scale structure (e.g., Cooray & Hu 2001b). 

Current predictions on statistics related to large scale 
structure weak lensing and their ability to measure cosmo- 
logical parameters are based on several assumptions; for 
example, it is implicitly assumed that background galax- 
ies, from which galaxy shape measurements are made, are 
uniformly distributed. Another assumption is the use of so- 
called Born approximation where one integrates along un- 
perturbed photon geodesies instead of perturbed photon 
paths. Here, we discuss the former assumption while Cooray 
& Hu (2002) presents a discussion of the corrections result- 
ing from dropping the Born approximation and including 
the so-called lens-lens coupling between two lenses at two 
different redshifts. 

The effect of background source clustering was first dis- 
cussed with respect to the lensing three point statistics, 
such as convergence skewness (Bernardeau 1998). Hamana 
et al. (2001) includes an extended discussion of this contri- 
bution. Recently, the same effect was revisited by Schneider 



© 0000 RAS 



2 Cooray 

et al. (2002) as a possible source of curl-like modes in lens- 
ing statistics. Here, we present a general discussion of the 
corrections resulting from clustering of background sources 
and consider effects two, three and four-point correlation 
functions. We suggest that the contributions are negligible 
and are unlikely to affect the current and upcoming weak 
lensing observations as a probe of cosmological parameters 
by comparing to predictions that neglect background source 
clustering. 



2 CALCULATIONAL METHOD 

We discuss our statistics in terms of lensing convergence, 



K{h) = -n m Hij dxrisix) I dx'g(x, x)<Hx'n; x) , (1) 



where n a (x) is the normalized radial distribution of back- 
ground sources such that f dxn B (x) = 1 an d <?(x'iX) is the 
lensing weight function when a background source is at a 
radial distance of x- 



(2) 



Here, x ls the radial distance, or lookback time, from the 
observer, given by 



dz' 



o H(z')' 

and the analogous angular diameter distance 



d A {x) = #o 



l n K 1/2 S mh(H n]( 2 x) . 



(3) 



(4) 



with the expansion rate for adiabatic CDM cosmological 
models with a cosmological constant given by 

H 2 = Hi [fi m (l + zf + Q K (1 + zf + Qa] . (5) 

Here, Ho can be written as the inverse Hubble distance to- 
day cHq 1 = 2997.9/i~ 1 Mpc. We follow the conventions that 
in units of the critical density 3H 2 /StvG, the contribution 
of each component is denoted Qi, i — c for the CDM, b for 
the baryons, A for the cosmological constant. We also define 
the auxiliary quantities Q m — fl c + fib andQn = 1 — ^\ Qi, 
which represent the matter density and the contribution of 
spatial curvature to the expansion rate respectively. Note 
that as Q,k — * 0, c(a(x) ~* X an d we define x( z = °°) = Xo- 
Though we present a general derivation of second order con- 
tributions to weak lensing due to background source clus- 
tering, we show results for the currently favorable ACDM 
cosmology with Q b — 0.05, Q. m = 0.35, Qa = 0.65 and 
h = 0.65. 

We assume that background sources are clustered in 
radial space such that n s (x) ~ ^s(x)[l + ^ n s (x)] where fluc- 
tuations in the number counts of background sources are 
related to that of the density field via a time, and possibly 
a scale dependent, bias 



5ns(x) = b 1 (x)S(x) + h 2 ( X )S 2 (x). 



(6) 



Here we have considered, perturbatively, contributions up 
to the second order in density perturbations. It is impor- 
tant that we consider contributions up to the S 2 term since 
these terms can contribute at the same second order level in 
the power spectrum. Though we discuss the effect in terms 
of lensing convergence, it should be noted that our calcula- 
tion equally well apply for, correlations of shear, such as the 
gradient, or electric- like, modes (see, Schneider et al. 2002). 

In Fourier space, we can write the first, second and third 
order contribution to convergence as 



K (1) (l) = -Q m H'l I dfBe 



"7 

Jo 

x £dx'g(x',x) J ^5(k, x V k ' x ' 9 , 



dx^six) 



« (2) (1) = ^ 



M J 



d 2 de~ ve 



dxn s (x)bi(x) 



and 



i0-(k lx +k 2X ') 



«< 3 >(i) = ^n, 



(7) 



,H 2 J d 2 0e- 19 £ 



dxn a (x)b2(x) 



d 3 ki 



<5(ki,x) 



(2tt)3 v ' a; 7 (2tt) 3 

f rf 3 k 3 /N-iS-Ckix+kax+ksx') 

J (2vf { A,X > 



5(k 2 ,x) / dx'g{x',x) 
o 



(8) 

respectively. 

We define the associated angular power spectrum, bis- 
pectrum and trispectrum of weak lensing convergence, in 
flat-sky as appropriate for current and upcoming experi- 
ments as 

(k(I)k(I')) = (2^) 2 M1 + 1')CT 
(k(1)k(1>(1")) = (27r) 2 fo(l + ... + l")B K (l,l',l") 
<k(1)k(1')k(1")k(1'"))c = (2^) 2 fc(l + ... + l"')T K (l, l', l", l'") 



(9) 



Here, kj...j = kj + . . . + kj and 8d is the delta function 
not to be confused with the density perturbation. Note that 
the subscript c denotes the connected piece, i.e., the trispec- 
trum is defined to be identically zero for a Gaussian field. 
Here and throughout, wc occasionally suppress the redshift 
dependence where no confusion will arise. 

We define the power spectrum of density fluctuations 



<<5(k)<5(k')) = (2nY i S(k + -k')P ss (k), 
where 



k 3 P(k) 



n+3 



T 2 (k), 



(10) 



(11) 



2tt 2 " \H 0/ 

in linear perturbation theory. We use the fitting formulae of 
Eisenstein & Hu (1999) in evaluating the transfer function 
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T(k) for CDM models. Here, 8h is the amplitude of present- 
day density fluctuations at the Hubble scale; with n = 1, 
we adopt the COBE normalization for 8h (Bunn & White 
1997) of 4.2 x 10 -5 , consistent with galaxy cluster abundance 
(Viana & Liddle 1999), with as = 0.86. To capture the non- 
linear aspects of the power spectrum, we use the prescription 
by Peacock & Dodds (1996). 

We will now discuss contributions to the angular power 
spectrum, bispectrum and trispectrum of convergence from 
clustering of background sources. To illustrate results, we 
take a redshift distribution for the background sources of 
the form 

M«) = d2) 

where (a, (5) denote the slope of the distribution at low and 
high z's, respectively around a mean- like parameter given by 
~ zq. For the purpose of this calculation we take a = j3 — 1.5 
and take zq to be 1.0 so as to mimic the expected back- 
ground sources from current and upcoming lensing catalogs. 
In general, since clustering evolves to low redshifts, with a 
decrease in the background source redshift distribution to a 
lower redshift, we expect an increase in the importance of 
second order effects associated with clustering of background 
sources. 

2.1 Power Spectrum 

To the first order, using (re' 1 - 1 (I)k^(I')), we simplify with the 
Limber approximation (Limber 1954) in Fourier space fol- 
lowing Kaiser (1992; Kaiser 1998) to obtain the well known 
result that 

I 2 



cr 



9„2 u 4 



* rf ../ g 2 (x'.x) p 



(13) 



The second order convergence power spectrum result- 



ing from source clustering involves two terms: (k 



(2)^(2) 



and (k^k^). The latter includes an additional, and equal, 
term through a permutation. We simplify these contribu- 
tions again using the Limber approximation (Limber 1954) 
such that 

( K (2) (l) K (2) (l')> = (27r) 2 fo(l + l') 



<?h p f h 



(2tt)2 \d A (x)' 



dA(x') 



(14) 



and 



(k (3) (1) k (1) (1')} = (27r) 2 fc(l + l') 

I 



d 2 h_ p ( h 



(2nr \d A ( X ) 



X 



dA( X ') 
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(15) 

Note that the total contribution to the convergence power 
spectrum follows as (k (2) (l)/t (2) (l')>+ (k (3) (1)k; (1) (l')}± 
(k (1) (1)k (3) (1')>. We denote the first contribution by C ; 22 
term and the latter two terms by Cf 1 . Thus, the total con- 
tribution to the power spectrum due to source clustering 
is 



C, 



cf 2 + cf 1 



or . ¥u*jy*S0L£«^ 



d 2 A (x' 



d 2 h 

(2±P 



Ii 



n 31 _ 9 2 rr4 



d 2 h 



d A ( x )' A J' \d A (x>) 



d 2 A (x) Jo d A ( X ') 



(27T) 



II 



d A (xY 



x)p 



I 



d A (x') 



(16) 

In figure |l|, we show the second order correction to the 
angular power spectrum of convergence. The solid line is 
the well known first order result, while the dashed line is the 
contribution from the Cf 2 term. Note that C 22 oc b 2 , and we 
have taken the value of bi = 1 for illustration purposes. The 
long-dashed line is the Cf 1 contribution, where Cf 1 oc 62 
and we have taken 62 = 1 for simplicity. Since we do not 
have detailed information on the galaxy bias, to illustrate 
our results, we have taken the bias to be redshift and scale 
independent. Adding a redshift dependent bias of the from 
b(z) ex (1 + z) 1 , however, did not lead to a significantly 
different result from the one suggested in figure [j] when | — y | < 
2. 

There is one aspect of bias that should be kept in mind 
when interpreting figure [l]. In general, quadratic bias is ex- 
pected to be negative, such that the two terms, Cf 2 and 
Cf 1 , added together will give a contribution which is lower 
than what one would naively expect if simply added to- 
gether. In the IRAS PSCz catalog, Feldman et al. (2001) 
finds I/61 = 1.20±g;i| and b 2 /b\ = -0.42 ± 0.19. Similar 
results our available from the 2dF survey by Verde et al. 
(2002): 61 = 1.04 ±0.11 and b 2 = -0.054 ±0.08. Since there 
is no conclusive evidence for a non-zero value for b 2 , the first 
approximation that b 2 — and 61 = 1 leads to the dashed 
line with the conclusion that second order effects are gener- 
ally below the few percent level for multipole less than 10 4 
corresponding to angular scales less than few arcminutes. 

On the other hand, if b 2 = —0.5 and 61 = 1, we obtain 
the dot-dashed line as the total contribution to the angular 
power spectrum of convergence due to source clustering; the 
shape of the power spectrum is due to the fine cancellation 
of Cf 2 and Cf 1 terms. We suggest that, in addition to the 
linear bias, the extent to which background source cluster- 
ing affects statistics such as power spectra or correlations 
depends on the detail aspects of galaxy biasing such as the 
quadratic bias. In any case, we find that source clustering 
effects are unlikely to be a strong contaminat for current 
lensing experiments. 
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Figure 1. The angular power spectrum of convergence. The solid 
line shows the first order contribution while corrections due to 
background source clustering is shown with dashed {Cf 2 ) and 
long-dashed (C ; 31 ) lines. We have assumed b± = 1 and b 2 = 1 
in these two cases, respectively. The dot-dashed line is the total 
second order power spectrum when b\ = 1 and b 2 = —0.5 con- 
sistent with suggestions in the literature for galaxy bias (see text 
for details). The dotted line is the C 31 contribution when galaxies 
trace the linear density field instead of the fully non-linear power 
spectrum. 



There is also another important aspect related to Cf 1 . 
The integral over li denotes the power spectrum traced by 
galaxies and Cf 1 effectively scales with this integral as an 
overall normalization. If galaxies do not fully trace the non- 
linear power spectrum, as predicted by the Peacock & Dodds 
(1996) formulae for the dark matter, then the contribution 
would be lower than what we have predicted. We can bracket 
the expected range of variation by replacing the power spec- 
trum involved with the integral over li with that of the linear 
power spectrum, as it is generally expected that the galaxy 
power spectrum lies between the linear and non-linear cases 
of dark matter. Since the contribution to the integral here 
comes from all angular scales, the bahvior of either the lin- 
ear or non-linear power spectrum at small angular scales 
becomes to some extent important for the calculation pre- 
sented here. As we do not have a reliable method to predict 
the non-linear power spectrum at small scales, we safely cut 
off the calculation at k ~ 10 6 h Mpc -1 . 



2.2 Bispectrum 

We can write the resuling contribution to the bispectrum by 
considering terms such as (k ( ' (1i)k (1) (12)k (1) (13)) and add 
the necessary permutations. We write one of these terms as 



,(2) 



(1)K 



(1) 



(1> (1) (!")> = (2tv) 2 S d (1 + 1' + 1") (^a m H- 



Jo 



4(x') 



dA{x'V 



Jo Jo 



d\(x') \dA( X ') ,X 



(17) 



so that the bispectrum is 

B K (h,h,h) = (^n m H^y 



rxo rx 

x / dxns{x) I d X 
Jo Jo 



d\{x') \d A (x'Y 



77, X +Perm. , 



(18) 



where the permutations are with respect to the ordering of 
(h,h,l3) and involves five additional terms. 

Following Cooray & Hu (2001a), we can write the third 
moment of convergence using the bispectrum as 

( K 3 (ff) ) = J_ ./'(2Ji + l)(2i2 + l)(2J 3 + l) 



47T 



h h I3 




B? lhl3 W h (a)W h (a)W h (a).. 



(19) 



where the quantity within parantheses is the Wigner-3j sym- 
bol, which in the case of no angular dependence can be writ- 
ten as 



h h h 




(-1) 



L/2 



(#)! 



(#-'i)K#-fa)K#-is)! 

1/2 



(L - 2h)l(L - 2l 2 )\(L - 2l 3 )l 



(L + iy. 



(20) 



for even L; it vanishes for odd L. We refer the reader to 
Cooray & Hu (2000) for additional details on the Wigner 
3-j symbol. 

Using the third moment, we construct the skewness as 



S3 (a) 



(21) 



where the all-sky expression for bispectrum, in terms of the 
flat-sky derivation, is 



/(2Ji + l)(2J2 + l)(2! s + l) / h h h 
- V 4^ I 



xB K (h,l 2 ,h) . 
Similarly, the second moment is defined as 

1 



(22) 



(23) 



where Wi are the multipole moments (or Fourier transform 
in a flat-sky approximation) of the window. For simplicity, 
we will choose a window which is a two-dimensional top hat 
in real space with a window function in multipole space of 
Wi(cr) — 2Ji(x)/x with x = la. 
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Figure 2. The skewness of convergence. We show the results 
based on gravitational perturbation theory calculation to the sec- 
ond order (PT; dashed line), and the calculation valid for highly 
non-linear regime based on hyper-extended perturbation theory 
(HEPT; Scoccimarro & Frieman 1999; dotted line). The solid line 
is the expected skewness based on the halo model calculation fol- 
lowing Cooray & Hu (2001). For comparison, we show results from 
N-body particle-mesh simulations by White & Hu (1999). The 
higher order correction to skewness resulting from background 
source clustering is shown with a dot-dashed line. The second or- 
der contribution is at the level of few tens of percent of the total 
expected from non-linear evolution of gravitational perturbations. 

In figure we summarize our results on the expected 
contribution to the convergence skewness. We refer the 
reader to Bernardeau (1998) and Hamana et al. (2002) for 
an extended discussion on the effects on background source 
clustering on skewness. As shown in figure ^| and discussed 
in prior publications, the contribution to skewness is at the 
level of few tens of percent and depends strongly on param- 
eters such as the mean redshift of background sources and 
the width of the redshift distribution. A higher mean red- 
shift and a smaller width result in a smaller contribution to 
skewness while a lower mean redshift and a broader distri- 
bution can contribute up to 30% or more of the skewness 
expected from non-linear evolution of gravitational pertur- 
bations. 
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Figure 3. The ratio of non-Gaussian to Gaussian error. The solid 
line is the calculation for non-linear evolution of gravitational 
perturbations using dark matter halos following the prescription 
in Cooray & Hu (2001b), while the dotted line is contribution 
arising from source clustering effects. We use the same binning 
scheme in multipole space as Cooray & Hu (2001b). 
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so that the trispectrum is 
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2.3 Trispectrum 

We can write the resuling contribution to the trispectrum 
by considering terms such as (re^ 2 - 1 (Ii)k^ (la)/*^ 1 ^ (Is)/?' 1 -* (U)) 
and add the necessary permutations. We write one of these 
terms as 

( K (2) (1) K (2) (1') K (1) (1")^ (1) (!'")> = 



(27r) 2 fc(l- 



Hr 



dx 



[bins] 2 



(Xi 



d 2 A (xi 



JI3+I1 

dA(xi, 



Xi)+P 



f |1 3 + b 

\d A ( X i) 



where the permutations are with respect to the ordering of 
(h,h,h, U)- 

For the purpose of this calculation, we assume that up- 
coming weak lensing convergence power spectrum will mea- 
sure binned logarithmic band powers at several h's in mul- 
tipole space with bins of thickness 5k. 



C 



dj_j_ 



«««(-!) 



(26) 



here A B (U) = J d 2 l is the area of 2D shell in multipole and 
1 itc.n'be written as A s (li) = 2nli5k + n(5li) 2 . 
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We can now write the signal covariance matrix as 



1 

A 



A s i 



T"CU, -li.lrf.-lj): 



(27) 
(28) 



where A is the area of the survey in steradians. Again the 
first term is the Gaussian contribution to the sample vari- 
ance and the second the non-Gaussian contribution. A re- 
alistic survey will also have shot noise variance due to the 
finite number of source galaxies in the survey. For a com- 
parison of previous calculations, we take the same binning 
scheme as the one used in Cooray & Hu (200fb) and used 
in White & Hu (1999). 

In figure 3, we show the ratio of -R = Ca/C^^ where 
C« aus is the contribution with simply the Gaussian variance. 
This ratio can also be written as 

AsiTa 



R= 1 + 



(29) 



and we plot R — 1 to highlight the difference between source 
clustering and non-Gaussian aspect of large scale structure. 
As shown, the clustering only leads to a few percent con- 
tribution, at I ~ 10 3 , beyond the Gaussian variance while 
the non-Gaussianities due to large scale structure clustering 
contributes at the level of 10% or more. One can safely ig- 
nore the relative increase in the variance of power spectrum 
measurements due to background source clustering. 

An additional aspect of the covariance resulting from 
non-Gaussianities is that band power estimates are corre- 
lated. These correlations can be written as 



a 



(30) 



In figure |i|, we show the behavior of the correlation coeffi- 
cient between a fixed lj as a function of three Ws. When 
li = lj the coefficient is 1 by definition. Due to the presence 
of the dominant Gaussian contribution at U — lj, the co- 
efficient has an apparent discontinuity between U — lj and 
li — lj-i that decreases as lj increases and non-Gaussian 
effects dominate. As shown, however, the correlation coef- 
ficients due to the non-Gaussian nature of the large scale 
structure is over an order of magnitude larger than than 
the correlations resulting from the clustering of background 
sources. The results related to the covariance suggests that 
non- Gaussian effects resulting from the clustering of back- 
ground sources is not expected to strongly influence the abil- 
ities of weak lensing experiment to obtain precision measure- 
ments of cosmology. 



3 SUMMARY & CONCLUSIONS 

We have discussed the second order contributions to weak 
gravitational lensing convergence resulting from the cluster- 
ing of background sources from which galaxy shape mea- 
surements are made in weak lensing experiments. The clus- 
tering of source galaxies induce a second order contribution 
to the two point statistics such as weak lensing convergence 
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Figure 4. The correlation coefficient Cij, as a function of the 
multipole li with lj, as shown in the figure. The top lines are 
the correlations due to the non-linear evolution of gravitational 
perturbations, as calculated under the halo approach of Cooray 
& Hu (2001b). The bottom lines are the resulting correlations 
due to the clustering of background sources from which shape 
measurements are made. These correlations are at the level of few 
percent, while non-Gaussian nature of the large scale structure 
induces correlations at the tens of percent level or more. 



angular power spectrum. For the angular scales of interest, 
we have shown that this contribution is at the level of few 
percent; to some extent, however, the exact contribution is 
uncertain due to unknown aspects associated with galaxy 
biasing such as the quadractic bias. 

Our calculations related to the skewness generated by 
background source clustering is consistent with previous cal- 
culations by Bernardeau (1998) and Hamana et al. (2001). 
We have discussed a new non-Gaussian aspect of the back- 
ground source clustering involving a contribution to the four 
point correlation function of shear or the trispectrum in 
Fourier space. The non-Gaussian four-point function is of 
interest since it determines the covariance of power spec- 
trum measurements. The background source clustering in- 
creases the Gaussian covariance at the level of few percent 
when I ~ 10 3 . This increase, however, is an order of mag- 
nitude or more below the increase resulting from the intrin- 
sic non-Gaussian nature of the large scale structure due to 
the non-linear evolution of gravitational perturbations. The 
trispectrum contribution to the covariance also leads to cor- 
relations between band power estimates, though, these are 
again at the few percent level or below and are unlikely to be 
a significant source of error for current and upcoming weak 
lensing experiments. 
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